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NUMERICAL ANALYSIS OF THE TRANSIENT RESPONSE 


OF AN AXISYMMETRIC ABLATIVE CHAR LAYER 
CONSIDERING INTERNAL FLOW EFFECTS 

By Claud M. Pittman and Lona M. Howser 
Langley Research Center 

SUMMARY 

The differential equations governing the transient response of the char layer of an 
ablating axisymmetric body, internal pyrolysis gas flow effects being considered, have 
been derived. These equations have been expanded into finite-difference form and pro- 
gramed for numerical solution on a digital computer. Numerical results compare favor- 
ably with simplified exact solutions. 

The complete numerical analysis was used to obtain solutions for two representative 
body shapes subjected to a typical entry heating environment. Pronounced effects of the 
lateral flow of pyrolysis gases on the mass flow field within the char layer and the asso- 
ciated surface and pyrolysis interface recession rates are shown. 

INTRODUCTION 

One-dimensional ablation analyses have been used extensively to study the thermal 
response of charring ablator heat shields subjected to aerodynamic heating. Recent 
experimental results, both in ground and flight test, show that a two-dimensional effect, 
involving the lateral flow of pyrolysis gases within the char layer, can significantly affect 
the performance of ablative materials. (See refs. 1 and 2.) This phenomenon becomes 
important in heat-shield areas such as stagnation regions and leading edges which have 
relatively sharp radii of curvature, high heating rates, and large surface-pressure gradi- 
ents. The high heating rates produce thick char layers and, because of the large surface- 
pressure gradients, the pyrolysis gases within the char layer tend to develop a significant 
lateral flow component. Because of this lateral flow, the mass of pyrolysis gas injected 
into the boundary layer is reduced in some areas and increased in others and this condi- 
tion affects both blocking of convective heating and surface recession. 

One -dimensional ablation analyses, such as that described in reference 3, can be 
modified with factors obtained from exact solutions for simplified two-dimensional flow 
fields. Such a procedure (which was used in ref. 1) is useful but does not provide any 



assessment of the accuracy of the results that are obtained. Clearly, a two-dimensional 
analysis of heat and mass transfer within a char layer is needed. Several two-dimensional 
thermal analyses are available. (See refs. 4 to 10.) Of these, only reference 9 attempts 
to treat the pyrolysis gas flow field within a charring ablator. In reference 9, only the 
energy and continuity equations are solved and the simplifying assumption is made that 
the gas mass flow rate vector in the char is oriented with the temperature-gradient vec- 
tor rather than with the pressure -gradient vector. 

This paper presents an analysis which determines both the temperature and pressure 
fields. within the char layer. The analysis considers (1) the effect of two-dimensional 
mass flow on the temperature distribution and recession rates and (2) the effect of tem- 
perature on the pressure and, therefore, on the mass flow distribution. The analysis is 
transient and is programed in terms of general curvilinear coordinates. 

The governing differential equations are derived and are converted to finite- 
difference form by using Taylor's series expansions. An acceptable method of solu- 
tion for the finite -difference equations was difficult to find. The method which was 
finally used to obtain solutions was to solve implicitly through the thickness of the 
char layer and explicitly along the other coordinate. This method is stable when the 
gradients through the thickness are dominant and yielded several solutions of interest. 

A method of solution which is always stable was not found. 

As the numerical solution is time consuming at best, and unstable under certain 
conditions, a criterion (based on a simplified two-dimensional exact solution and a one- 
dimensional numerical analysis) is presented which can be used to identify conditions 
where lateral internal mass transfer will be significant. 

Numerical results compare favorably with exact solutions for simplified cases. 
Solutions are also shown for two body shapes, subjected to a typical ballistic entry envi- 
ronment, to illustrate the significant effects of two-dimensional mass flow. 

SYMBOLS 

A quantity defined by equation (3) 

Ai,A s constants in mass loss rate equations corresponding to specific reaction rate 
B dummy variable used in appendix D 

Bj,B s constant in exponential of mass loss rate equation corresponding to activation 
energy 
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c 


constant used in exact solution (see eq. (Dl)) 


C e oxygen concentration by mass at edge of boundary layer 

Cp specific heat of char 

Cp specific heat of pyrolysis gases 

f porosity of char 

h enthalpy 

^l’h2’hs scale factors given by equations (A5) 

Ah c heat of combustion 

Ahp heat of pyrolysis 

K reaction- rate constant for char oxidation 

k thermal conductivity of char 

L total distance along base curve (fig. 2) 

M molecular weight of gases 

m mass flow rate 

p pressure 

q heating rate 

R radius of curvature 

R 0 universal gas constant 

R^,R 2 radii (defined after eq. (D8)) 

R s normal distance from axis to base curve (fig. 2) 



radius (defined after eq. (D8)) 


x,y ,e 


V 

K 

X 


p 

p 

p' 

a 


dimensionless radius (eq. (D8)) 

arc length 

temperature 

velocity 

coordinates 

distance from base curved to pyrolysis interface 

absorptance of char surface 

weighting factor for transpiration effectiveness 

angle between tangent to base curve and axis of symmetry (fig. 2) 

angle defined in figure 3 

char thickness 

emittance of char surface 

transformed y- coordinate 

permeability 

mass of char removed per unit mass of oxygen 

viscosity 

density of char 

density of uncharred material 
Stefan- Boltzmann constant 
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T 


time 


cp quantity defined by equation (12) 

angle in spherical coordinates (eq. (D8)) 

Subscripts: 

C cold-wall convective 

C,net hot-wall convective heating corrected for transpiration 

c char 

e at edge of boundary layer 

g gas 

o initial 

* 

p pyrolysis interface 

R radiative 

s char surface 

t stagnation point 

w wall 

x x- component 

y y-component 

1- D one-dimensional 

2- D two-dimensional 

An arrow over a symbol denotes a vector. 
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ANALYSIS 


A transient analysis is developed for an axisymmetric char layer (fig. 1) subjected 
to a typical entry heating environment. The primary result sought is the effect of axi- 
symmetric temperature and pressure distributions on the mass flow within the char layer 
and on the char surface and pyrolysis interface recessions. Although only the char layer 
was analyzed, the back surface of the uncharred material was used as the base curve for 
the coordinate system, since it always remained fixed. 

Assumptions 

The primary assumptions used in the analysis are as follows: 

(1) The density of the char layer is constant. 

(2) The thermal conductivity of the char layer is a function of temperature only. 

(3) The permeability of the char layer is a function of temperature and direction. 

(4) The heat transfer and pyrolysis gas flow are zero normal to the surfaces x = 0 
and x = L. (See fig. 2.) 

(5) Char surface removal is by oxidation only. 

(6) The effect of surface shape changes on surface heating rates and pressures 
are neglected. 

(7) The gradients in the 0-direction are zero. 

Differential Equations for the Temperature Distribution 

Governing temperature equation .- The governing equation for the temperature dis- 
tribution in the char layer in vector form is 



Heat Heat absorbed by Heat 

conducted pyrolysis gases stored 


The development of equation (1) into the form to be used in the analysis is given in 
appendix A. Through this development, equation (1) becomes 
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i_/^2 k 9 t\ _ d_(H kA 8T\ _ A k 32^ + A -i^kA ^ 

9x J ax^hj 817 J 9T?\hl 9 xj drjj 


+ ly( h l h2k If') " h2lh x fe( c p T ) + h 2 Am x ^( c p T ) 




where 


. . st h l h 2 c p^c s T h 1 h 2 (l - ^mpPCp aT 

" h l h 2 pC P9 T + 6 9 tj + 6(p’-p) 9r] 


( 2 ) 





(3) 


y - y P 
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(4) 


and and h 2 are scale factors, shown in figure 3 and discussed in appendix A. 

Initial conditions .- The initial conditions that must be specified are the temperature 
distribution, mass-transfer rates, and the body shape. For most cases of interest, the 
initial temperature distribution is uniform and the initial mass-transfer rate is zero. 

Boundary conditions .- The boundary conditions on equation (2) are as follows: 

At 77 — 0 (the pyrolysis interface), the heat conducted in is equal to the heat 
absorbed by the pyrolysis reaction 


k 8 T 
8 drj 



Ah p 


where rhp is obtained from an Arrhenius -type equation 
rhp = Aje T 


(5) 


The effect of the uncharred material on the ablative performance of the char layer can be 
roughly approximated by including the sensible heat absorption of this material in the heat 
of pyrolysis Ahp. 

At 77 = 1 (the char surface), the heat conducted in is equal to the net heat input to 
the surface 
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( 6 ) 


k 8T 
6 8r? 


= % 


ere T 4 


where (see ref. 3) 




1- 0.6 ^ s * y + 0.084^ as m y j 


+ aq^ + m c Ah c 


The quadratic expression in my is a second-order approximation to a curve obtained 
from the boundary -layer program of reference 11 for air to air injection. The expres- 
sion is valid for both blowing and suction at the wall (positive and negative my) . The 
mass loss rate of char at the surface m c is obtained from an equation for first-order 
oxidation (see ref. 3) and is 


where 


and 


^ _ K Pw C e 

IIl£ — 

j KPwC^e ~ hw) 
Xq C,net 

.Es 

K = A s e T 


1 - 0.6 ^a s m y + 0.084^ as m. 



(7) 


At x = 0, and x = L, the temperature gradient is zero 


DT _ 9T » 9T _ o 
dx dx~ 9rj 


( 8 ) 


Although at x = L, this boundary condition is usually physically incorrect, it should be 
adequate especially if this boundary occurs outside the area where two-dimensional 
effects are significant. 


Differential Equations for the Pressure Distribution 

Governing pressure equation .- To obtain the governing differential equation for the 
pressure distribution in the char layer, it is assumed that the pressure distribution can 
be described by using Darcy's Law (ref. 12) 


v = 



(9) 
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Multiplying both sides of equation (9) by p and assuming p„. = yields 

& & IVq X 

ih = - ^ Vp 2 


2pR Q T 

Then by conservation of mass, 


( 10 ) 


-V • m = V • cpV p 2 = f 


Dp 


9t 


(ID 


where 


= m 

V 2pR 0 T 


( 12 ) 


and f is the porosity of the char layer. 

By using the same development which was used for the temperature differential 
equation, the pressure equation (eq. (11)) becomes 



= hi h 2 f 


8pg 

9t 


me 8p g 

6 p Sri + 


(1 - Tj)m p 9p g 
(p’ - p)6 *T 


(13) 


Initial conditions .- The initial conditions that must be specified are the pressure 
distribution and the body shape. The initial temperature distribution and mass loss rates 
that are required to solve the pressure equation (see eqs. (12) and (13)), are obtained from 
the solution of the temperature equation. 

Boundary conditions .- The boundary conditions on equation (13) are as follows: 

At rj = 0 (the pyrolysis interface), the pressure gradient is proportional to the 
pyrolysis rate 

(Hi 

dri cp 
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At rj = 1 (the char surface) 


P 2 (x) = P e 2 (x) 

where p e is specified as a function of location and time. 

At x = 0, and x = L, the mass flux across the boundary is zero 


Dp 2 8p2 Sp 2 

Prv ^ m X " 


3x 9x “ dri 


( 15 ) 


(16) 


The previous discussion of the temperature boundary conditions at x = L applies also 
to the pressure boundary conditions at this location. 

Mass flow equations .- When the pressure distribution in the char layer has been 
determined from equations (13) to (16), the flow field m x and my can be determined 
from equation (10). Thus, 


m x = - 



- A 


8pj 
3 1 ] 


m y = "7r 


c o 3p 2 
6 3 7 ] 


(17) 


Coordinate Singularity 

Equations (2) and (13) apply to the entire region of interest; however, a line of dis- 
continuity (singularity) exists. An inspection of equation (A 5) shows that at x = 0, that 
is, along the line of symmetry, the scale factor hg vanishes. This coordinate singu- 
larity can be eliminated by using proper approximations valid only near x = 0. The 
derivation of the differential equations which are used at the singularity is given in 
appendix B. The temperature equation at x = 0 is 



9 T hiCpr?m c grj, hj&pPCp(l - rj) gip 

l pc P 9r + 6 3r? + 5(p ! . p ) 8 tj 


(18) 
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The pressure equation at x = 0 is 


9 JlIsL 9 8 (v k 8 P 2 ^ _ 9 A J_/iL 

4 ax^hj ax j ' 4 ax y*! 977 y " A 877 \h x ax 




ap2\ <p 9 p% 
877 + WE 877 


9 p g hjfrjihc 9 p g hif(l - 77 )m p 9 p g 
1 dr 5 p 9 r? 6(p’ - p) 977 


( 19 ) 


NUMERICAL PROCEDURES 


Finite-Difference Equations 

The differential equations were converted to finite-difference form through the use 
of the Taylor's series expansions. Forward, central, and backward differences were 
used. A summary of the finite-difference equations which are used in this analysis is 
given in appendix C . * 


Method of Solution 

It was extremely difficult to find an acceptable method of solution for the finite - 
difference equations. Both standard matrix inversions and explicit methods were tried 
but were abandoned when the stable time interval became so small that 1 hour of compu- 
puter time produced no significant progress. 

The equations could not be solved entirely implicitly because of the nonlinear terms 
and the interdependence of the equations. Thus, an iterative procedure was required. The 
first choice of an implicit method was the alternating- direction method. (See ref. 13 .) 

This method was not successful because, although the solution in the 77-direction (columns) 
was reasonably accurate, the solution in the x-direction (rows) was not. It is postulated 
that the reason for this result is that the zero- slope boundary conditions at x = 0 and 
x = L (eqs . (8) and ( 16 )) do not serve to bound the solution. When only derivative boundary 
conditions are used, the shape of the curve can be determined, but the magnitude at any 
point cannot. 

Other implicit methods of solution which were tried included a modified alternating- 
direction method (ref. 14 ) and a diagonal alternating-direction method (ref. 15 ). The 
method which was finally used to obtain solutions was to solve only the column equations 
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implicitly. The x-derivatives appear explicitly in the column equations. This approach 
is successful as long as the 77- derivatives are considerably larger than the x-derivatives. 
When the lateral flow m x becomes very large, some instability in the solution results. 

If the calculation continues to run, the instability generally subsides and does not appear 
to contribute significant errors in the final solution. In more severe cases, the solution 
diverges sufficiently to cause the program to stop or to make the solution meaningless. 

A method of solution which always remains stable was not found. 

RESULTS AND DISCUSSION 

The numerical analysis was compared with three simplified exact solutions to deter- 
mine the accuracy of the finite -difference expressions. The exact solutions involve either 
the temperature equations or the pressure equations separately. The derivations of the 
exact solutions are given in appendix D. The complete numerical analysis was used to 
provide solutions for two body shapes subjected to a typical entry heating environment. 
Pronounced effects of the lateral flow of gases on the mass flow field within the char layer 
and the associated surface and pyrolysis interface recession rates are shown. 

Comparisons With Exact Solutions 

Three comparisons were made between the numerical analysis and exact solutions. 
The temperature equations are compared with an exact Solution for a rectangular char 
layer undergoing quasi- steady- state ablation. The pressure equations are compared with 
two steady- state exact solutions, one for a rectangular char layer and one for a hemisphe- 
rical section. 

Quasi-steady-state temperature solution .- The comparison between the quasi- steady- 
temperature distribution obtained with the exact solution (eq. (D3)), and the distribution 
from the numerical analysis, is shown in table I. The input values used are shown, and 
the percent error at representative locations is given. The percent error values are cal- 
culated from the equation shown in the table. The maximum errors occur at the x = 0 
locations and decrease as x increases. The maximum error of 0.277 percent is con- 
sidered to be acceptable. The average error is 0.162 percent. 

Pressure distribution for a rectangle .- The comparison between the pressure dis- 
tribution for a rectangle obtained with the exact solution (eq. (D6)) and the distribution 
from the numerical solution is shown in table II. The input values used are shown and 
the percent error at various locations is given. The maximum errors occur at large 
values of x and small ?j- values. The largest error is 0.18 percent and the average 
error is 0.065 percent. 

The errors in the mass -flow- rate components are shown in table HI. The errors 
are an order of magnitude larger than the errors in the pressure distribution. The mass- 
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flow-rate components (eqs. (D7)) are obtained from derivatives (eqs. (17)), and it is very 
difficult to obtain accurate derivatives numerically. It appears that these errors are 
within the accuracy of the numerical method. Efforts to decrease the errors by including 
higher order terms in the derivatives were unsuccessful. 

Spherical pressure distribution .- The comparison between the pressure distribution 
for a hemisphere obtained with the exact solution (eq. (DIO)) and the distribution from the 
numerical solution is shown in table IV. The input values used are shown and the percent 
error at various locations is given. The largest error is 0.794 percent and the average 
error is 0.289 percent. 

The errors in the mass -flow- rate components (eq. (Dll)) are shown in table V. The 
errors are again an order of magnitude larger than the pressure errors. 

Typical Calculated Results 

Typical calculated results are given for two entry body shapes; these results were 
obtained by using the enthalpy, stagnation heating rate, and stagnation pressure histories 
shown in figure 4. These histories are typical of a ballistic entry from earth orbit. 

The two body shapes used in the calculations are shown in figure 5. Figure 5(a) 
shows a hemispherically blunted cone shape and figure 5(b) shows a hemisphere -cone 
body. The part of each shape considered in the analysis lies between the x = 0 and 
x = L lines. Beyond the x = L line, the surface conditions become nearly constant 
with x. The material property values used in the calculations are given in table VI. 

Hemispherically blunted cone .- The heating rate and pressure distributions around 
the blunt-faced body are shown in figure 6 as a function of distance along the reference 
surface. (See fig. 5(a).) Figure 6 shows a sharp pressure gradient at the first tangent 
point. (This point is the point where the large nose radius intersects the small corner 
radius.) Because of the large pressure drop around the corner, a low-pressure area 
exists in the vicinity of a high-pressure area. This condition can produce a large lateral 
flow effect because the gases within the char tend to flow laterally toward the low-pres- 
sure region rather than in the direction normal to the pyrolysis interface. 

The effect of the lateral flow within the char layer is shown in figures 7 and 8. Fig- 
ure 7 shows the mass flow, normal to the surface, at a point immediately preceding the 
sharp pressure drop. This station, in the analytical solution, is where the largest lateral 
flow effect is observed. A calculation with lateral flow and one in which lateral flow is 
not allowed (m x = 0) are shown. In the early part of the trajectory, no lateral-flow effect 
is evident because the lateral-flow effect is dependent not only on the pressure gradient 
but also on the value of the pressure and on the char thickness. These latter two values 
are both small early in the trajectory. 
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As the char thickness and pressure increase later in the trajectory, lateral flow 
not only decreases the mass flow at the surface but actually produces an inflow of the 
boundary- layer gases into the char layer. The decrease in mass flow out through the 
surface decreases blocking of the convective heating. When the surface mass flow 
becomes negative, the convective heating rate becomes higher than it would be for an 
impermeable surface. 

Figure 8 shows total surface and interface recession as a function of distance along 
the reference surface. Both the lateral-flow solution and the no- lateral-flow solution 
are shown. There is a very small difference in recession at the stagnation point x = 0. 
The lateral-flow effect is inversely proportional to the nose radius. Thus, although the 
effect is present, it is small because of the large hemispherical radius of the blunted 
cone. The differences in recession near the corner are significant. There is more 
recession at the corner than at the stagnation point for both cases because of the heating- 
rate peak. (See fig. 6.) However, the case with lateral flow shows considerably more 
recession. The increased surface recession is caused by increased oxidation which is 
due to the higher convective heating associated with reduced blocking and eventual inflow 
at the surface. The more significant increase in pyrolysis interface recession is caused 
mainly by the increased temperature gradients caused by the decrease in char cooling by 
the pyrolysis gases and the eventual influx of hot boundary- layer gas into the char layer. 

Near x = L, the lateral flow case shows less recession than the no- lateral flow 
case because the gases that were flowing laterally through the char come out through the 
surface in this region of low pressure, increase the blocking of convective heating, and 
decrease the surface temperature. 

The high recession rate near the corner has been observed experimentally on abla- 
tive heat shields of this shape. Previously, this large recession has been attributed to 
mechanical removal. The results presented here show that the lateral flow effect with 
its accompanying high oxidation rates can cause the observed shape change. 

Hemisphere-cone .- The heating rate and pressure distributions on the hemisphere- 
cone body (fig. 5(b)) are shown in figure 9. Figure 10 shows the mass flow normal to the 
surface at the stagnation point (x = 0). When the char- layer thickness and pressure 
become large enough, the mass flow through the surface decreases and eventually becomes 
negative. This condition causes the increased recession shown in figure 11. Again the 
difference in recession is greater at the pyrolysis interface because the influx of hot 
boundary- layer gas has a greater effect on the internal temperature gradient than the 
increase in convective heating has on oxidation. 
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CRITERION FOR ASSESSING THE IMPORTANCE OF LATERAL FLOW 


Since the analysis described in ibis paper is difficult to use, it is desirable to have 
a criterion which will indicate whether lateral flow should be considered or whether a 
one-dimensional ablation analysis is adequate. This criterion can be developed from the 
exact solution for the pressure distribution in a hemispherical char layer and a one- 
dimensional ablation analysis such as the analysis of reference 3, A lateral-flow param- 
eter for a hemisphere can be obtained from appendix D. From appendix D (eq. (D12)), 
the mass flow at the surface is 


R 


\p,r 


Ri 


m v = 


R 

-3«V(cos2^ - 


R 


^,r 


r 2 

Rl 


2 + af5i' 3 

3\R 2 , 


( 20 ) 


The point on a hemisphere where lateral flow has the largest detrimental effect is at the 
stagnation point. At this point xj/ = 0° and R^ r = R 2 . By assuming that the char thick- 
ness is small so that terms of order 52 can be neglected, equation (20) can be reduced 
to 


rily 



( 21 ) 


If this mass flow rate is small compared with the pyrolysis rate of a one -dimensional 
stagnation-point analysis, then the lateral-flow effect should not be significant. 

To determine at what value the mass flow of equation (21) can be considered small 
compared with the pyrolysis rate in a one- dimensional analysis, the calculations for the 
hemisphere-cone (fig. 5(b)) were evaluated. At 40 seconds the ratio of the mass flow of 
equation (21) to the one-dimensional pyrolysis rate is -0.001, at 50 seconds the ratio is 
-0i012, at 60 seconds the ratio is -0.145, and at 70 seconds the ratio is -0.65. An exami- 
nation of figure 10 shows that the lateral-flow effect becomes important between 50 and 
60 seconds. Therefore, it appears that if 


( m y>p2-P 

(H-D 


< 0.1 


a one-dimensional analysis is probably sufficient. 


( 22 ) 
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Although this criterion was developed for a hemisphere with constant char thickness, 
it should be useful for most blunt bodies if the char thickness is small compared with the 
radius of curvature of the surface. For example, if the corner radius of the blunted- cone 
shape (fig. 5(a)) is used in equation (21), then equation (22) predicts that the lateral-flow 
effect becomes significant after 50 seconds. This result is consistent with the results 
shown in figure 7. 

CONCLUDING REMARKS 

Differential equations governing the temperature and pressure fields within an axi- 
symmetric char layer were derived for assessing the effects of two-dimensional mass 
flow within the char layer. The equations were written in terms of general curvilinear 
coordinates, expanded into finite -difference form, and programed for numerical solution 
on a digital computer. 

Numerical results compare favorably with exact solutions for simplified cases. 
Numerical solutions for the complete set of equations are shown for two body shapes sub- 
jected to a typical ballistic entry environment. Results show that including the effect of 
the lateral flow of gases within the char layer has a significant effect on the performance 
of the ablation material. In areas where the lateral flow effect is pronounced, blocking 
of the convective heating is reduced and the mass flow normal to the surface may even 
become negative. In these areas there is a significant increase in the surface tempera- 
ture and in the surface and pyrolysis interface recession. The results shown through 
this analysis may explain some experimental ablation material recession data which 
could previously only be attributed to mechanical removal. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 28, 1972. 
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APPENDIX A 


DEVELOPMENT OF GOVERNING DIFFERENTIAL EQUATION 

This appendix describes the derivation of the differential equation that will be con- 
verted into finite-difference form. 

The governing vector equation for the temperature distribution in the char layer is 

7 ■ ( : M = If^p + 'v r p) T ] (A1) 

Heat absorbed by Heat 

pyrolysis gases stored 

By assuming that the density of the char layer is constant, that the specific heats of 
both the char and the pyrolysis gases are functions of temperature only, and that Pg « p, 
the right-hand side of equation (Al) can be written as 

If[( pc p + p g'p) T ] * ^P !£ + ! P T T* 

Equation (Al) can now be written as 


V • kVT - 


Heat 

conducted 


V • kVT - m • v(c p t) = pc p p- + c p T (-^ + v * *) 

Since the quantity within the parentheses in the last term is the continuity equation, 
the last term is equal to zero. Therefore equation (Al) becomes 

V • kVT - m • v(c p Tj = pCp p (A2) 

Because the surface and the pyrolysis interface locations are changing with time, 
it is convenient to make the transformation 

y - y p 

73 = ~ 5 ~ 
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APPENDIX A - Continued 


where 

y p thickness of uncharred material 

6 thickness of char layer 

With this transformation, the 77 -coordinate always goes from 0 to 1 through the char thick- 
ness 5. The result is a moving coordinate system in the y-direction. 

Because 6 and y p are functions of x, 77 is also. Therefore, total derivatives 
rather than partial derivatives must be used in the x-direction. This total derivative can 
be written (ref. 16) as 

D( ) a JL 

9x 9x 9x 877 “ 9x " drj 

where 

A 9x 6 v 7 9x 9 X / 

Total derivatives with respect to time must also be used to give 

DT = 9T 9 n 9T 

9t 9 T 9t 9tj 

and since 

9tl_9t 7 96 877 8 yp 

9t “ 96 9t + 9y p 9 t 

then 


The change in 6 is 

6 . 6o = -^dT + J^dT 


DT _ 9T , 

9r 9 t I 96 


91 + 

9r 9y. 


8 ?7 8 yp\ 


9T 

dr J drj 


(A3) 
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APPENDIX A - Continued 


where (see fig. 2) 

m c mass loss rate due to char removal at surface 

m p mass flow rate at pyrolysis interface due to pyrolysis reactions 

From the preceding equation, 

36 _ A c . m P 
3t P p' - p 

The change in yp is 



Thus 

8 Yp _ A p 
8t p' - P 


Then, since 

£2? _ V 
36 6 

877 1 

8 Yp " 6 

equation (A3) becomes 

DT 3T 
dr 3r 

The equation is transformed to a curvilinear coordinate system through the use of 
scale factors. (See ref. 10.) The differential arc length ds in the curvilinear coordi- 
nate system shown in figure 3 is 

(ds)2 = h^^(dx)2 + h 2 ^(d<p )2 + hg^dy)^ 


V m c (1 - J))m p 

6P (p' - p)6 


3T 

9 T) 


(A4) 
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APPENDIX A - Concluded 


where the scale factors are 


hl = 1 + L 


h 2 = R s + y cos j3 / 


h 3 = l 


J 


(A5) 


The curvilinear coordinate system should conveniently describe any axisymmetric 
body geometry of interest. However, if the Cartesian coordinate system is required, 
care must be taken to make all scale factors equal one. 

To make the analysis more tractable, it is assumed that no temperature gradients 
exist in the 0-direction. This assumption limits the analysis to axisymmetric shapes and 
flow fields. 

Through the use of the preceding development, the final form of the differential 
temperature equation is obtained by combining equations (A2), (A4), and (A5) and is 


_ 9 _ 

8x 







+ A 




+ 



fe( c P T ) + h 2 Am x J?( C P T ) - 


h l h 2 rh y a 
6 3?J 



= h h , pc 8T + h l h 2 C P T ' m c |T ^ h l h 2» - qHPCp 9T (A6) 
1 2 p 8t 6 8 v (p- _ p)8 8 n 

The procedure used to develop the pressure equation is very similar and will not 
be given. The final form of the differential-pressure equation is given by equation (13). 
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APPENDIX B 


DERIVATION OF THE DIFFERENTIAL EQUATIONS 
AT THE COORDINATE SINGULARITY 

In this appendix, only the derivation of the temperature equation at x = 0 will be 
given. The derivation of the pressure equation at the coordinate singularity is similar. 

The derivation of the temperature equation when h 2 = 0 is initiated by determining 
the behavior of h 2 as x — 0. As x — 0, R g — x and cos ~ Therefore, 

( h 4-0 ■ 4 + 1) <BI 


and 


l d M 

\ 9x / x-0 



(B2) 



(B3) 


The temperature equation is evaluated at the coordinate singularity by using equa- 
tion (A2) with the scale factors included. Equation (A2) can be written as 



The first two terms of equation (B4) must be modified because of the coordinate sin- 
gularity. By using equations (Bl), (B2), and (B3) and L' Hospital’s rule, the first term can 
be evaluated as follows: 


1 _9_/^2 k 9T] = J_ JL/JL 8T.\ + k 9T 9h 2 _ 1 _8_/k_ 9T\ k 9T 
h l h 2 9x l h i 9x ) h i 9x lhj 9x ) hj2h 2 9x 9x h l 9x (hi 9x ) h 1 h 2 9x 
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APPENDIX B - Concluded 


Then, since |^-= 0 at x = 0, 

7 Civ 7 


_9_/k_ 92^ 

k 8T lim 9x l h l 8x / iWlE) 

h l h 2 9 X x— 0 9h 2 /9x h l 8x \ h l **) 


Therefore, 


1 8 ft 

^1^2 ^ \h 


2 k |T 

i 9x 


'-■x— 0 


2_ _8_/k_ 9T’ 
h^ 9xlh^ 8x 


(B5) 


A similar procedure is used to evaluate the second term 


i-JL/ hlhok |T 




= JL -§_/hik 2X\ + — 2. = 1- e_4 li k + k 8T x 

h x 9y \ ! 9y j + h 2 9y 8y h A 9y f 1 9y ) + + y_j 8y R 


Therefore, 


1 8 


hjh 2 9y 


8T 

9y 


x— 0 


= _L ±. h , k m ,±_dT 

hi 9y\ 1 9yJ + h 1 R9y 


(B6) 


The differential equation for the temperature distribution at the coordinate singu- 
larity can now be given. By using equations (B5) and (B6) in equation (B4) and including 
the transformation of equation (3), the temperature equation at x = 0 is 


9 9 /k 9T\ o 8 ZkA ST^ 2 A 9 /k 8^ , 9A 9 /kA 8T\ 
2— (——I - g -J ~ 2A ^ gjj + 2A J 


Sxlhj 0x 


_8 _a 

llk ST^ 

+ x. 

9T _ •„ 

iz. 9_/c t\ 


1 97/) 

R6 

877 

1 6 9?7 ( C P / 


h i pc P ft + 


8T h l c ?V*Q 9T ^(l “ ^“P PC P 8T 


Br] 


(p f - p) 6 977 


(B7) 
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APPENDIX C 


FINITE -DIFFERENCE EQUATIONS 
Temperature Equations 

The finite-difference temperature equations are as follows: (The m subscript 
corresponds to the 77 - coordinate, m = 1 being rj = 0 and m = s being rj = 1, The 
n subscript corresponds to the x coordinate, n = 1 being x = 0 and n = L being 
x = L.) For 1 < m < s, 1 < n < L, the equation is 
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APPENDIX C - Continued 


For m = 1, 1 < n < L, the equation is 


•k T 


■ A,» 4‘- n+1 


( ho \ /ho \ / ho \ 

M-k + Uk T in+ r k T 1 n 1 

i hi / , \ h-. / i A > n \h-, / A >n-i 

\ 1 /l,n+l \ 1 /l,n-I \ 1 /l,n-i 


57 A6rb P Ah P , - (eT A6rh P Ah p) t 
k 1 A.n+i ' 1 


57 k ) ( T 2 ,,rVl) (kN ( T 3 ,n + 1 - T 3 ,n-l) 3 ( ET k ) ( T l,n + 1 " T l, 

/ 2 ,n ' \ 1 / 3 , 1 a \ 1 /l,n 


A r) Ax 


4 At? Ax 


4 A?7 Ax 




3 Ar;^ 


07h7 A6rh P Ah P 


( T 2,» - T l,») 


( T 3,n - T 2.n) 


* k 3 ( h 'H „ ’ A,* - ( h * h2k )3, n 3A,2 - 


(c p T) - (c T) 4 (c ■ T) - fc T) - 3 /c T) 

(^ 4 ,„ 2,n 2 aT ( j ‘-“ 

( h l h 2 m y) 1 , n ( CpT ) 2>n ( P )s,n 3 ^ pT \,n _ [ hih 2 PCpmp1 4T 2,n ' T 3,n ~ 3T l,n 
6 n 2 At 7 " (p' - p) 6 , 2 A?7 


( h l h 2 pc p) 1 n 


^ - ( T o) 1>n 
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APPENDIX C - Continued 


For m = s,l<n<L, the equation is 
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APPENDIX C - Continued 


For 1 < m < s, n = 1 , the equation is 


6 3 ( T m,2- T m,l) 5 ( T m,3- T m,2) 


3 (Kl) 3 ( T m+l,l ' + ( T m+1,2 ' T m-1.2 


2 At? Ax 


i hi ) 5 ( T m+1,2 ' ^m-1,2) + ( T m+1,3 T m-1.,3 > 

V m, 2 


6 Ax At? 


t J( hlk )m + Il T m + l,l '( h l k )m + l,l + ( h » k )m-Il Tm > 1 . Mm-l/”- 1 ' 1 


Or ^ L 


T T My) (W -(c p T) 

T m+l,l • T m-l,l 1 ^m,!^ Wl,l \ V-1,1 


R ’%i,l 2AT? 5 1 


h l c ^c\ ^m+l,! ~ (1 ~ ? ?)^l” 1 pP c p ^m+l,! ~ ^m-1,1 


(p’ - P)6 j 2 AT] 

J m.l 




' ( T °)m,] 
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APPENDIX C - Continued 


For l<m<s, n=L, the equation is 



3 (ir kA ') 

w u.i 


T m+1,L “ T hi-1,l) + ( T m+1,L-1 ” T m-l,L-l) 


4 A t] Ax 


H 3 

* 1 /m,L-| 


' T m+1,L-1 “ T m-1,L- 1) + (' 


12 A?7 Ax 


T m+l,L-2 “ T m-l,L-2) 





At/** 


2A?7 


fa h 2 c p 77m c \ T m+ l |L ~ T m-1,L 
' 6 m,L 2 Ar ? 


^ 1-77 m p pc P h l h 2 j T m+1,L " T m-1,L 
vP’- p 5 / msL 2 Af? 


“ ( h l h 2P c p) 


■’m,L " (T 0 ) 


m,L 


m,L 


At 


= 0 


(C5) 
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APPENDIX C - Continued 
For m = 1, n = 1, the equation is 



Ax 
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APPENDIX C - Continued 


For m = s, n = 1, the equation is 
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APPENDIX C - Continued 
For m = 1, n = L, the equation is 



_ '^hgpcpihp j 4T 2,L ~ T 3,L ~ 3T 1,L 
(P* - P)5 , 2 Ar l 

x T 1.L ' ( t °)i,l „ 

- (“l^j L At 0 (c8) 
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APPENDIX C - Continued 


For m = s, n = L, the equation is 



3 Ax 2 




6 L (q t «T 4 ) +5 L _i(q t - « t4 ). 


+ 6 t 

's,L L ' 
2 Ax 


s,L-l 





6 Ax 


+ 



( h l h 2 k ) S -|,L( T s-l,L - T s-2,l) 


3 Aiy* 



+ 0?( h 

J 

( h 1 h 2 lh y) S ,L ( g P T )s-2,L * 3 fi T )s,L ~ 4 ( r P T )s-i,L 
6 l 2 At? 


/hiVpm^) 3T s,L + T s- 2,L ~ 4T s-l,L 




's,L 


2 At? 


(h lh2 pc p ) Sj] 


t s,L ' ( T o) SjJ 


At 


= 0 


(C9) 


31 



APPENDIX C - Continued 
Pressure Equations 

The finite-difference equations for the pressure distribution in the char layer are 
as follows: 

For l<m<s, 1 < n < L, the equation is 
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APPENDIX C - Continued 


For m = 1, l<n<L, the equation is; 
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APPENDIX C - Continued 


For m = s, 1 < n < L, p 2 = p2 , For 1 < m < s, n = 1, the equation is 

s,n e S) n 



Ax^ 3 Ax^ 



,2 

m+1,1 


- P 


2 

m-1,1 


+ P 


- P 2 

m+.l,2 m-1,2 


2 Arj Ax 



hl W“ 


p2 - p2 \ + fp2 _ p2 

m+1,2 m-l,2y l m+1,3 m-1,3 


6 Ax Arj 


( h 1 ^). m+i 1 P m+ 1 , 1 


1 + 
m+-,i 




P m,l + ( h l^)m-i 


^l 2 


Arj* 


<p \ P m+l,l P m- 1,1 /hjf? 7 m c ^ ( Pg )m+l,l ( P ®)m-l,l 


l R 6/ m 1 2 ^ 


6p 


/m,l 2 


hjf(l - T])m p 

6(p' - p) 


-*m,l 


Pg Wl,l " \ S/m- 1,1 
2 A?? 


M 


( Pg )m,l " ( P g,o) m ,l 


m,l 


At 


= 0 


(C12) 
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APPENDIX C - Continued 


For 1 < m < s, n = L, the equation is 




- p 2 


+ p: 


-p2 


m+l,L m-l,L/ I m+l,L-l m-l,L-l 


4 A ?7 Ax 


1 /m,L-| 


2 2 

m+l,L-l ^m-ljL-li rm+l,L-2 r m-l,L-2 


+ P 


12 At? Ax 


( h l h 2 <p ) m+ l L P m +l,L 
1 2 


, hlh 2%4L + ( h l 


m-|,L 


P 


m,L 


*L* 


At ?' 


At?' 


( hl Hm-i,L P m-l,L 
At ? 2 


?m r 


/hjhgfT?] 

\ 6p / T 

' 'm,L 


Per) - T "" (Pg) 1 T 
m+l,L \ & /m-l,L 

2 At? 


h 1 h 2 f(l - 7?)m p 


6(p' - p) 


m,L 


f ^ Li,L " ri-i, L 

2 At? 


(w) 


( P e)m,L ■ (Mm.l 


m,L 


At 


= 0 


(C13) 
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APPENDIX C - Continued 


For m = 1, n = 1, the equation is 





(C14) 
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APPENDIX C - Continued 
For m = 1, n = L, the equation is 



h l h 2 flfa P 4 ( Pg )2,L " ( Pg )3,L ~ 3 ( Pg )l,L 
6(p’ -p)j 2 At? 

i,Jj 




(Cl 5) 
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APPENDIX C - Continued 


Mass- Flow Equations 

The finite-difference equations for the mass-flow distribution in the char layer are 
as follows: 

Component in x-direction.- For 1 < m < s, 1 < n < L, the equation is 


Nm,n 


f (p\ 

P 2 i " P 2 i 

m,n+l m,n-l 

A m,n( p m+l,n " p m-l,n) 

(hj 

2 Ax 

2 Ar) 

\ Vm,n 

■ 



(C16) 


For m = 1, I < n < L, the equation is 



+ A 


l,n 



For m = s, 1 < n < L, the equation is 


(C17) 



For 1 S m = s , n = 1, and L, the equation is 

(»x) , = (™x) ~ 0 

\ /m,l \ /m,L 

Component in y -direction.- For 1 < m < s, 1 = n = L, the equation is 


(Cl 9) 



a P “ - p" i . 

y m,n m+l,n m-l,n 

6n 2 Ar] 


(C20) 
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APPENDIX C — Concluded 


For m = 1, 1 ln = L, the equation is 


H,a = N, 


For m = s, 1 § m < l, the equation is 


m 3/p 2\ _ 4p2 + p2 

\ _ ^s,n \ e /n s-l,n s-2,n 

l m 4,n = " "5n 2 Tr, 


(C21) 


(C22) 
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APPENDIX D 


EXACT SOLUTIONS 

Quasi-Steady-State Temperature Solution 

An exact solution is obtained for a rectangular char layer undergoing quasi-steady- 
state ablation. Steady-state environmental conditions and constant material properties are 
used. Conditions are selected so that the char thickness is constant both in the x-direction 
and with time. Mass flow in the y-direction only is considered. 

Under these conditions, the governing differential equation for the temperature dis- 
tribution is 


a 2 T a 2 T c 9T 
8x2 8y 2 9y 


0 


where 


T 


= T-T 0 


(Dl) 


C 



Constant 


The boundary conditions used are 



T = 0 



cos 


2 flx 

2L " 2 \ 1 


+ cos — 
L 


The exact solution is 


(x = 0; x = L) 


(y = 0) ) (D2) 


(y = 5) 


_Cy 

q e ^ 

H o 

^sinh SL 

... 2 

cos ^ sinh ^ ) 

_C6 
ke 2 

_C6 
\ce 2 

B cosh ^ _ c sinh 


(D3) 
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APPENDIX D - Continued 


where 



Pressure Solution for a Rectangle 

An exact solution is obtained for the pressure distribution in a rectangular char 
layer with a steady-state surface pressure distribution and constant temperature. Under 
these conditions, the governing differential equation for the pressure distribution is 


8 2 p 2 8 2 p 2 
— — + — = 0 

3x“ 3y ■ 

The boundary conditions used are 



The exact solution is 


(D4) 


(x = 0; x = L) 

(y = 0) > (D5) 


(y = 6) J 



+ cos 


cosh 


iry\ 

L 


cosh 


/ 


(D6) 


The mass flow components are also calculated from the exact solution. Thus, 


m x = 


„ 2 p 2 cosh 

8x 2 L L cosh £» 

Li 


m v = -<p 


8p Z 

9y 


p, 2 

_ (D _L L cos 2^ 
^ _ cos _ 


sinh 
cosh ^ 


(D7) 


41 



APPENDIX D - Concluded 


Spherical Pressure Distribution 

An exact solution is obtained for the pressure distribution in a constant-thickness 
char layer on a spherical quadrant. A steady-state pressure distribution and constant 
temperature are used. Under these conditions, the governing differential equation for 
the pressure distribution is 


_9_ 

9r 


(r2 — ^-|sin P 

\ 9 r / sin ip dxp\ 


aq-o 


r 2 


where r is a dimensionless radius R^^/Rj and 1 g r I g® 


The boundary conditions used are 


(D8) 


9pr 

dip 

dp 2 

W 


= o 


0 


p2 = p^cosfy 
The exact solution is 




0; rp 
(r 

It- 


= 1) > (D9) 


R2 

R 


1 1 


p2 = ^ + P t 2 (c°s2^-|^- 

Ri 


+ 2 ( R 1 'f 


3 R 




+ 2pl 

3\R 2 y 


(DIO) 


The mass flow components are calculated from 


m x " m p ~ " R 


(p 9p 2 dcp o . , \ R 1 

— -rj- = r- P/srn ^ cos xp ±— 
xp,r \ >r 1 


+ 2/%J 3 

3\R, 


W/ 


R af + 2/^A 3 

R l/ 3 r 2/ 


(Dll) 


(R. 


p,r 


Ri 


my = m r s= 


~<P 


9p^ 


9R. 


j//,r 




R' 




r 2 ' 

R iy 


+ 2P 

3 R 2 y 


(D12) 
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TABLE I.- COMPARISON OF EXACT AND CALCULATED 
QUASI- STEADY- STATE TEMPERATURE DISTRIBUTIONS 

[input values used were q = 113,5 kW/m2; k = 0.624 W/m-K; 

6 = 0.06096 m; L = 0.3048 m; C = 3.2214 1/m; T Q = 556 K; 
Arj = 0.25 (m = 5); Ax = 0.0381 m (n = 9)] 


Percent error = — — ^ x 100 

T exact 


V 


Percent error for 

x of ■- 


0 

0.25 

0.5 

0:75 

1.0 

1.0 

0.242 

0.237 

0.215 

0.143 

0.009 

.75 

.275 

.264 

.228 

.129 

.023 

.5 

.274 

.259 

.209 

.104 

.020 

.25 | 

.277 

.251 

.183 

.075 

.025 

! 

0 

.230 

.200 

.120 

.050 

-.020 




TABLE II.- COMPARISON OF EXACT AND CALCULATED 
PRESSURE DISTRIBUTIONS FOR A RECTANGLE 

[input values used were p g = 105.336 kN/m2; 6 = 0.06096 m; 
L = 0.3048 m; Arj = 0.25 (m = 5); Ax = 0.0381 m (n = 9^j 


Percent error = ^ ^ ca ^ c x 100 

y exact 


V ; 


Percent error for 

x’ of - 


0 

0.25 

0.5 

0.75 

1.0 

1.0 


---- 

----- 



.75 

0.000 

0.015 

0.000 

-0.076 

-0.070 

.5 

.009 

.025 

.000 

-.125 

-.180 

.25 

.014 

.061 

.000 

-.142 

1 

OO 

0 

.048 

.036 

.000 

1 

-.140 

-.172 




TABLE EL- COMPARISON BETWEEN EXACT AND CALCULATED 


MASS- FLOW-RATE DISTRIBUTIONS FOR A RECTANGLE 
\cp = 7.76 x 10" 11 m3/N-s] 


Percent error = ^exact ->calc x 100 

m exact 

(a) Percent error; x- component 


V 


Percent error for 

x of - 


0 

0.25 

0.5 

0.75 

1.0 

1.0 

- 

2.63 

2.69 

2.64 

-- 

.75 

- 

2.49 

2.49 

2.50 

-- 

.5 

- 

2.44 

2.41 

2.45 

-- 

.25 


2.42 

2.36 

2.42 


0 

| 

2.42 

j 

2.35 | 

i 

2.42’ 

-- 


(b) Percent error; y-component 

Percent error for x of — 
0.5 
0.000 
.000 
.000 





TABLE IV.- COMPARISON BETWEEN EXACT AND CALCULATED 
PRESSURE DISTRIBUTIONS FOR A HEMISPHERE 

{input values used were p g = 47.88 N/m.2; = 0.06096 m; R2 = 0.3048 m; 

6 = 0.24384 m; L = 0.0957 m; At] - 0.25 (m = 5); Ax = 0.0119 m (n = 9)] 


Percent error = exa ^ x 100 

exact 


V 


Percent error for 

x of - 


0 

0.079 

0.157 

0.236 

0.314 

1.0 

.75 






0.059 

-0.14 

-0.115 

-0.029 

-0.024 

.5 

-.638 

-.619 

-.352 

.174 

.477 

.25 

-.794 

-.738 

-.390 

.086 

.287 

0 

.028 

1 

-.045 

-.154 

-.281 

-.359 
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TABLE V.- COMPARISON BETWEEN EXACT AND CALCULATED 
MASS FLOW DISTRIBUTIONS FOR A HEMISPHERE 


[<p = 7.76 x 10" 11 m3/N-s] 
^exact “ “calc 


Percent error = 


X 100 


m exact 

(a) Percent error; x- component 



Percent error for 

x of - 


0 

0.079 

0.157 

0.236 

0.314 


-0.418 

-.550 

-3.17 

-6.13 

2.64 


0.075 

-.359 

- 2.88 

-5.8 

.676 


0.568 

.441 

-2.30 

-4.92 

4.1 


(b) Percent error; y-component 

Percent error for x of - 


0 


-4.17 

3.06 

4.44 

-.932 


0.079 


-0.890 

2.66 

2.86 

-.676 


0.157 


-0.180 

3.32 

3.39 

.273 


0.236 


-2.03 

1.46 

1.80 

-1.46 


0.314 


-3.05 

1.87 

2.74 

-1.31 























TABLE VL- MATERIAL PROPERTY VALUES USED IN CALCULATIONS 
FOR BODY SHAPES SHOWN IN FIGURE 5 


k, W/m-K, for - 

278 K 

833 ....... 

1111 

1389 ...... 

1667 ...... 

1944 ...... 

2222 ...... 

2500 

2778 

3056 ...... 

Cp, KJ/kg-K, for - 
2 78 K ..... . 

833 

1111 ...... 

1389 ...... 

1667 ...... 

1944 ...... 

2222 . 

2778 . . . . . . 

3333 ...... 


0.156 

0.156 

0.499 

1.217 

1.872 

2.652 

3.744 

4.742 

6.240 

7.675 

. 3.64 
. 3.35 
. 8.37 
13.60 
. 7.53 
. 4.39 
. 5.02 
. 9.20 
20.08 


c p , KJ/kg-K 
Ah p , MJ/kg 
A s , kg/m^-s 
B s , K . . . 
A^, kg/m^-s 
B i5 K . . . 


X 

€ ..... 
Ah c , MJ/kg 
fi , N-s/m2 

M 

K y ,K x , m 2 . 
p, kg/m3 . 
p', kg/m 3 . 
f 


.... 2.26 
.... 2.32 

4.88 x 10 10 
4.25 x 10 4 

7.71 x 1-6® 
1.29 x 10 4 
.... 1.0 
. . . 0.232 
.... 0.75 
. . . . 0.8 
.... 10.0 
4.79 x 10” 5 
. . . . 20.0 
9.29 x 10” 12 

.... 256 
.... 552 
.... 0.8 
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Figure 2.- Coordinate system 
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Base curve 
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Figure 5.- Axisymmetric shapes analyzed. All dimensions are in centimeters. 










blunted cone. 



No lateral flow 



Figure 7.- Comparison of mass flow normal to surface, with and without lateral flow, at point 
of maximum two-dimensional effect. x/L = 0.78; blunted cone. 




Figure 8.- Maximum surface and interface recessions. Blunted cone 



Heating rate and pressure ratios 




Figure 10.- Comparison of mass flow normal to surface, with and without lateral flow, at point 
of maximum two-dimensional effect, x = 0; hemisphere-cone. 
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